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ABSTRACT 



Context. The knowledge of the regular (large scale) component ef the Galactic magnetic field gives important information about the 
structure and dynamics of the Milky Way, as well as constitutes a basic tool to determine cosmic rays trajectories. It can also provide 
clear windows where primordial magnetic fields could be detected. 

Aims. We want to obtain the regular (large scale) pattern of the magnetic field distribution of the Milky Way that better fits the 
polarized synchrotron emission as seen by the WMAP satellite in the 5 years data at 22 GHz. 

Methods. We have done a systematic study of a number of Galactic magnetic field models: axisymmetric (with and without radial 
dependence on the field strength), bisymmetric (with and without radial dependence), logarithmic spiral arms, concentric circular 
rings with reversals and bi-toroidal. We have explored the parameter space defining each of these models using a grid-based approach. 
In total, more than one million models are computed. The model selection is done using a Bayesian approach. For each model, the 
posterior distributions are obtained and marginalised over the unwanted parameters to obtain the marginal (one-parameter) probability 
distribution functions. 

Results. In general, axisymmetric models provide a better description of the halo component, although attending to their goodness- 
of-fit, the rest of the models cannot be rejected. In the case of disk component, the analysis is not very sensitive for obtaining the 
disk large scale structure, because of the effective available area (less than 8% of the whole map and less than 40% of the disk). 
Nevertheless, within a given family of models, the best-fit parameters are compatible with those found in the literature. 
Conclusions. The family of models that better describes the polarized synchrotron halo emission is the axisymmetric one, with 
magnetic spiral arms with a pitch angle of x 24°, and a strong vertical field of 1 fiG at z « 1 kpc. When a radial variation is fitted, 
models require fast variations. 

Key words. Magnetic Fields - Polarization - Galaxy: structure 



1. Introduction 

Spiral galaxies exhibit large-scale magnetic fields. The Milky 
Way is not an exception, but obtaining its spatial distribution is 
extremely difficult. Most methods for observing magnetic fields 
' are b ased either on ra dio observations of the synchrotron emis- 
; sion (IWolleben et all [2006; Reich, 2006; Testori et al., l2008l 
and references therein), or on the Faraday Ro t ation (here i nafter , 
_ FR) of pulsa r s (e.g . ,Weisberg et al.L l2004t iHan et all 120061: 
. iNoutsos et al.L l2008l) and extragalactic radio sources (here- 
. inafter, EGRS) (e.g. 'Gaensler et al.', '2007; 'Br own et al.L l2007t 
[Haverkorn et al., 2008; Carretti et al., 2008). In addition, sev- 
eral radio lines show Zeeman splitting, which can be used as 
well to directly constrain the str ength of the magnetic field (e.g. 
iFish et aLLlIoollHan & Zhangl l2007). 

Despite the large numbers of pulsars and EGRS for which 
the Rotation Measure (hereinafter, RM) has been recently deter- 
mined, there is no consensus about the large-scale pattern of the 
Galactic Magnetic Field (GMF). Probably, i t is more comple x 
than previously e xpected, a s pointed out in iMen et al ] (I2008h . 
Recent results by Sun et al.] (1200 8) show an axisymmetric disk 
distribution with reversals inside the solar circle as the best 



model to describe the GMF when all-sky maps at 1 .4 GHz from 
DRAG and Villa Elisa, the 22 GHz map from WMAP satellite, 
and the Effe lsberg RM survey o f EGRS are combined. Results 
derived from Bro wn et al.l (l2007h by using RM of EGRS suggest 
an axisymmetric pattern of the disk magnetic field. Vallee ( 200^ 
claimed for an inclusion of a ring model to describe the field. 
The interest on the structure of this large scale GMF is justified 
by a number of reasons. First of all, this field might be of impor- 
tance when considering the dynamics of the galaxy at those large 
scales (Nelson, 1988; Battaner et al., 1992; Battaner & Florid^, 
[19951 iKutschera & Jalochail2004tlBattaner & Floridd 120071) . A 
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good characterization of the large-scale GMF pattern would 
allow a detailed correction of the galactic contribution for 
a better understanding of cosmological magnetic fields (see 
Battaner & Florido, 2009, for a recent review), which would 
be potentially obs erved with upcoming CM B missions, such as 
PLANCK ( iThe Planck Collaborationll2006h . 

In addition, the GMF modifies the trajectory of 
high energy cosmic rays, therefore its knowledge is 
crucial to understand their distribution in energy and 
direction, such as obtained in experime nts as Auger 
( Bluemer & for the Pierre Auger CoUaborationi 2008), Milagro 
dAbdo et all l2009l) and others. The global anisotropics found 
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by MiLAGRO have been interpreted by iBattaner et alJ (l2009l) as 
produced by galactic magnetic fields. Moreover, GMF could be 
important to explain the well known knee in the spectrum of 
cosm ic rays energies around 10'' GeV dMasip & Mastromatteol 
12008 ). as sub-knee-energy cosmic protons trapped by GMF 
should have a much larger optical depth for interactions with 
WIMP'S. 

The extraction of structure of the galactic magnetic field 
from the measurements of the polarized intensity is extremely 
difficult. The galactic magnetic field cannot probably be de- 
scribed by a single model but there co uld be at least three com- 
ponents, the thin disk, the thick disk ( Beuerm ann et al.L Il985h 
and the halo, each one defined by its own model and parameters. 
Taking into account that each component is largely unknown, an 
analysis of a multi-component magnetic field renders an objec- 
tive estimation of the best configuration a too complicated task. 

The thin disk is characterized by the highest field strengths 
and we are embedded in it. However, magnetic fields in the 
thin disk are predominantly random, and therefore they barely 
contribute to the observed net polarized emission, even if a z- 
component of the regula r fields exist (Ha n et al., 1999). On the 
other hand, local spurs ([Berkhuiisen et al.l Il971b . as the North 
Polar Spur, highly distort the main field configuration. Even at 
this high frequency Faraday depolarization cannot be neglected 
in particular regions. Observing at high galactic latitudes the thin 
disk field is contaminated by the thick disk and the halo fields, 
which could even become dominant. Models for the thick disk 
are scarce in the literature but those proposed for the thin disk 
could be tested, even if different parameters could characterize 
both disks. Our results for the disk inferred from the polarized 
synchrotron emission will give a complementary insight on the 
field structure. The halo structure remains unknown and has a 
very different structure consisting in a double torus in two hemi- 
spheres with opposite directions. After s ome pioneer detections 
(ISimard-N ormandin & Ki-onberg. 1980; Han & Oiao, 1994), it 
has been modelled by Han et al. (1991), Harm et 9.1. (1999.) . 
iTinvakov & TkachevI (l2002h . iProuza & Smfdal (l2003h and oth- 
ers. To illustrate the current uncertainty on this component, we 
can compare the maximum value of the magnetic strength, 1 jjG 
following Prouza and Smida, and 10 jiG following Sun et al., 
although these last authors propose a maximum strength of only 
2 fxG when the thermal electron scale height is increased by a 
factor of 2. The contribution of the halo field to the polarized 
emission is therefore difficult to estimate. Most models take only 
into account RM from extragalactic sources to estimate the halo 
structure. 

In this work, we carry out a systematic comparison of a num- 
ber of possible GMF models, exploring which one is providing 
a better fit to the large-scale polarization map at 22 GHz. 

Our analysis is b ased on the 5-yea r WMAP data release, and 
extends the work bv'Page e t al.l (l2007b by considering a detailed 
comparison not only with the polarization angle, but with the 
polarized intensity (i.e. Stokes's Q and U parameters). Although 
the polarization angle can be used to described some properties 
of the large scale pattern of the GMF, it is not sensitive to some 
parameters (e.g. the field strength) and it also contains some in- 
trinsic degeneracy with respect to the direction of the field lines. 
Because of this reason, our main results will be obtained with 
the analysis of the {Q, U) maps, although the independent anal- 
ysis based on the position angle will be done in some cases for 
comparison. 

In this work, we have used eight disk models and one halo 
model. Different masks enable us, in an indirect way, to esti- 
mate the different contributions of the galactic components at 



different galactic latitudes, but the consideration of a multicom- 
ponent field has not been fully undertaken. In what follows, we 
will speak of the disk without specifying if thin or thick. 

For obtaining the emission we also need a model of the distri- 
bution and spectrum of cosmic rays which is another important 
source of uncertainties. Here we have assumed that the cosmic 
ray structure follows that of the gas, as they are produced by su- 
pernovae and cannot travel faraway from the birth place. This 
assumption is rather usual in the literature but rather question- 
able too. 

Here we focus on the study of the large scale field. In princi- 
ple, we could divide the galactic magnetic field into three com- 
ponents: a random component for scales lower than say 100 pc 
(see Haverkornet al., 2008), for which some works h ave pub- 
Ushed the turbulence spectrum ( Han et all 120041: iHanl 12009); a 
main "spiral" field for scales typical of spiral arms, and a "galac- 
tic scale" main field. Some a uthors cons ider that 1 kpc is a length 
defining the large scale (e.g. lHani 1200 8). therefore taking spiral 
arms as a large scale phenomenon. In fact, the field direction 
is opposite in arms and in inter-arms regions (e.g. IBeck et all 
119961 ; iHan et a n l2006h . Here, however, we are interested in 
scales of the galaxy itself, being therefore spiral arms consid- 
ered as wavy perturbations. It is clear that we need to investigate 
these three components, but this separation is important because 
the tools, and mainly the interpretation in terms of the gener- 
ation mechanism, could be completely different. The random 
field should be produced by turbulence in a magnetized medium, 
supernova explosions and other local mechanisms. Spiral arms 
produce characteristic motions that enhance magnetic fields in a 
high conductivity medium. Fields at the galactic scale would be 
interpreted in terms of galaxy formation and/or dynamo effects. 

The polarized synchrotron emission is a valuable tool for in- 
vestigating the overall field pattern at large scales and specially 
at high galactic latitudes. The total emission is much affected by 
random fields and by non-polarized galactic emission as free- 
free (see e.g. Miville- Deschenes et all 12008). FR of pulsars is 
a powerful technique but it is very much affected by enhance- 
ments and tangling of the field by the passage of spiral waves. 
FR of extragalactic sources would inform about the larger scale 
fields but here the poor knowledge of the intrinsic FR is a heavy 
problem. When using all-sky observations such as provided by 
WMAP, and PLANCK in the future, the observed emission is 
integrated along a path traversing the galaxy, therefore the wavy 
effect of spiral arms is, at least in part, smoothed. 

This could explain why RM of pulsars and WMAP po- 
larization have given somewhat inconsistent results. There are 
other all-sky surveys at lower frequencies (see iReicR l200 6l) but 
Faraday depolarization and dust emission are very strong at lat- 
itudes below 30 deg. The detailed study of the WMAP polariza- 
tion measurements in terms of the galactic scale magnetic fields 
is more a complementary than an additional technique. It is rel- 
atively free of Faraday depolarization, of dust contamination, of 
random fields and of spiral waves perturbations. 

To complement our study, for some of the models we also in- 
vestigate if a modification of the radial variation of the strength 
of the GMF is producing an impact on the quality of the fit. This 
radial variation is also largely unknown and may be of impor- 
tance on the production of the all-sky polarization maps. 

The paper is organized as follows. In Section 2 we describe 
the set of GMF models that have been investigated in this anal- 
ysis. In Section 3, we describe the numerical method used to 
produce the template maps, as well as the the model selection 
method used to determine our preferred model and the corre- 
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sponding set of parameters. We present the results in Section 4 
and discuss them in the following section. 



2. Models of Galactic Magnetic Field 

We have focused our analyses on the constraints derived from 
the polarized intensity map at 22 GHz. At these frequencies, the 
physical process which dominates the polarized intensity is the 
synchrotron radiation emitted by the population of relativistic 
cosmic r ay (CR) electrons w ith energies between 400 MeV and 
25 GeV (IStrong etal.Ll2007h . as they interact with the GMF. 

Therefore, in order to obtain a prediction of the polariza- 
tion pattern at this frequency, we first require a description of 
the distribution of the relativistic CR electrons in our Galaxy. 
For the purposes of this paper, in which we are interested only 
in the large-scale pattern of the GMF, it is sufficient to con- 
sider a simplified description of the CR electron population. 
H ere, we use the spatial dis tribution of relativistic electrons as 
in lDrimmel & Spergell(l2001i) . i.e. 



Ne - No exp 



5 kpc 



sech 



1 kpc 



(1) 



where A^o ~ 3.2 x 10^'^cm"-' is derived from the value for the 
CR electron density on Earth (Sun et al., 2008), and r and z are 
the radial and the vertical coordinates in cylindrical galactocen- 
tric coordinates, respectively. This model essentially assumes the 
same spatial density distribution for the CR electrons as for the 
interstellar gas. At first order, this is what we would expect, be- 
cause the higher the interstellar gas density, the higher the star 
formation rate, and therefore, the higher supernova production 
which gives a higher relativistic electron density. These cosmic 
electrons lose energy (via synchrotron itself) in short distances 
of less than 1 kpc. 

The value for A^o in equation ([T]l is very uncertain. This value 
is usually obtained by assuming that the CR electron spectrum 
can be described by a power-law with constant spectral index 
p - T>. However, observations in the last fe w ye ars, as well 
as numerical simulations (see IStrong et all |2007l and refer- 
ences therein) suggest that this assumption is not appropriate 
for the entire spectrum. Moreover, different observations show 
variations of the order of 50% (or even larger) for this quantity. 
Therefore, we expect that this uncertainty might introduce a bias 
in the recovered amplitude of the field strength for the different 
models, and will be taken into account as explained below. 

In the following sub-sections we present the set of GMF 
models we investigate in this work, all of them taken from the 
literature. Most of these models have been proposed/constrained 
using the analysis of Faraday Rotation of pulsars and EGRSs. 
Thus, it is also interesting to explore if these models are also able 
to describe the large-scale polarization pattern seen in WMAP 
22 GHz maps. In total, we investigate eight models describ- 
ing the disk and halo fields, and one model conceive to de- 
scribe the halo field. The disk models are: 1) Axisymmetric 
(ASS), 2) Axisymmetric with radial dependence of the strength 
(ASS(r)), 3) Bisymmetric positive (BSS+) , 4) Bisymmetric neg- 
ative (BSS_), 5) Bisymmetric positive with radial dependence 
of the strength (BSS+(r)), 6) Bisymmetric positive with radial 
dependence of the strength (BSS+(r)), 7) Concentric Circular 
Rings (CCR) and 8) Logarithmic Spiral Arms (LSA). Moreover, 
for the halo component, we shall consider the bi-toroidal (BT) 
model. Throughout this section, all coordinates (r, z) refer to 
cylindrical galactocentric coordinates. 



2.1. Axisymmetric Model 

The a xisymmetric model (see e.g. IValleel 119911; iPoezd etakl 
11993 1) is one of the simplest descriptions of the GMF. It is com- 
patible with a non-primordial origin of the galactic magnetism, 
based on the dynamo theory. There are several possible models 
of the family of the axisymmetric distribution. The components 
for this model of GMF are given by: 



Br = Bo(r) sm(p) cosOf (z)) 
B^ = Boir) cos(p) cos(x(z)) 
B, = Bo{r)sm(x(z)) 



(2a) 
(2b) 
(2c) 



where p is the pitch angleQ which it is considered constant; B()(r) 
is the field strength (which in principle might be a function of the 
radial distance), and xiz) corresponds to a "tilt angle". For our 
study, we have adopted the following functional dependence: 



X(z) =^otanh(— ) 
Zo 



(3) 



where we use a value of zo - 1 kpc for the characteristic scale of 
variation in the vertical direction. 

For the computations in this paper, we shall consider two 
cases for the radial dependence of Bo(r). The first case (here- 
after ASS) corresponds to B{){r) = Bq (constant value). As a 
second case (hereafter ASS(r)), we consider a radial variation of 
the type: 

Bo(r) = (4) 

where ri represents the characteristic scale at which Bo{r) de- 
creases to half of its value at the galactic centre. This radial varia- 
tion is based on a possible extension of the model by Poez d et al.l 
il992)- The expression has appropriate asymptotic behaviors, in 
the sense that we obtain a finite value when r is close to the 
galactic center (r — » 0), and asymptotically tends to oc 1/r when 
r — > oo, as suggested in Battaner & Florido (2007). We must 
note that in equation (|4|, Bi and ri are not independent, provided 
that we fix the value of the GMF strength in the solar neighbour- 
hood. For example, using Rq - 8 kpc for the Sun galactocentric 
distance, and Bq - 3 fiG for the magnetic field strength in the 
solar neighbourhood, we can re-write equation ^ in terms of a 
single free-parameter as: 



Bo(r) = 



3ri + 24 
ri + r 



(5) 



where r is given in kpc and Bo(r) in fiG. 

Summarizing, a particular ASS model is fully described once 
these three parameters are given: [B(),p,xq], while the ASS(r) 
model, one would in principle require four parameters. However, 
for the ASS(r) family of models, we will use the constraint given 
in equation (|5]l, which in practice means that we will only have 
three free parameters to represent a certain model: [ri,p,xo]- 

The typical range of variation of the p values found in the 
Uterature for the ASS model of the disk (see .Valleej, 1199 Ih 
is shown in Table I . In general, different pitch angle and field 
strength values are derived for the different spiral arms of the 
Galaxy. 



' The "pitch angle" is defined here as the angle between the azimuthal 
direction and the magnetic field. The azimuthal direction (0) increases 
in anti-clockwise, so p is positive since the anti-clockwise tangent to the 
spiral is outside the circle with radius r. Note that in some works in the 
literature, the convention for the ^-angle is exactly opposite (i.e. it in- 
creases clockwise). Finally, note that for the case p = 0°, the solenoidal 
model is recovered. 
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2.2. Bisymmetric Model 

This model is compatible with a primordial origin of the cosmic 
magnetism. It could explain the reversals of the magnetic field 
derive d from obse rvations o f RM of pulsars (e.g. Han & Oiao, 
Il994t iHanl 120011: iHan et al.L ^Mi . The three components for 
this model are given by: 



Br = BoWcos(0+^lnl — ))sin(p)cosO^(z)) (6a) 



= Boir) cos 1 + ;8 In I — 1 1 cos(p) cos(x(z)) (6b) 



B,^Bo(r)sm(x(z)) 



(6c) 



where Bo(r) is the field strength: = 1 / lan(p), being p the pitch 
angle: is the distance Sun-Galactic center (=s 8 kpc): and fi- 
nally xii) is the tilt angle, which we assume is also given by 
equation ([3]). Note that in Eq. (|6]l we are considering two pos- 
sible families of bisymmetric models, which correspond to the 
positive (negative) sign inside the big parenthesis for Br and 
B^. In this way, we are considering the two conventions that 
are found in the literature. "Positive" bisymmetric models (here- 
after, BSS+) are defined with the same convention to recover the 
model proposed by Han & Qiao ( 1994), while "negative" bisym- 
metric models (here after, BSS+) use the same sign convention of 
iJansson et al.l ('2009'). Note that for BSS_, the spiral has opposite 
sign for the magnetic direction. 

As in the previous case, we consider two sub-families of 
models. The first one, noted as BSS±, corresponds to constant 
strength of the magnetic field. In this case, a model is com- 
pletely defined by giving three parameters: [Bq, p,xo]- The sec- 
ond family, noted as BSS±(r), includes a radial variation of the 
field strength according to equation (|4]i. For this case, we again 
fix the field strength to the value in the Solar neighbourhood 
(equation|5]l, so a model is specified by giving three parameters: 
[ri,p,Xo]- 

Typical parameter values for the disk model found in the 
literature, after converted to our convention for the sign of the 
pitch angle, are p - (8.2 + 0.5)°, = 1.4 yuG and B global - 
(1.8 + 0.3) nG (IHan & OiaoL [1991: v = 14° and Bq = 1 yuG 
(ISimard-N ormandin & Kronber gi ll980l) : and p = 11° and B,, = 
(2.1 ± 0.3) nO (.Han etal... 2006,) . 

2.3. Concentric Circular Ring (CCR) model 

This model was proposed bv iRand & Kulkarnil (Il989h to fit the 
RM of their pulsars catalogue, with the aim of providing a way to 
take into account reversals of the magnetic field at diff'erent radii. 
In their original expressions, they do not consider a vertical com- 
ponent for the magnetic field (i.e. B, = 0). Here, we shall extend 
the eq uations for this model presented in llndrani & Deshpandd 
(119991) 
way: 



to account for a vertical dependence, in the following 



Bo 



Br^ 

^ sin(7rDr/w) 
= BosinOr(z)) 



sin ( 1 cosCr(z)) 



(7a) 
(7b) 
(7c) 



where w is the space between reversals: Dr is the distance to the 
first reversal: Bo is the field strength, and finally ;i^(z) is given by 
equation All distances (w and Dr) are given in kiloparsecs. 
Note that we introduce an additional factor sin(7rD, /w) in the 



definition of B^, so Bo still preserves the meaning of the mag- 
netic field strength in the solar neighbourhood. 

In our analyses, we have not considered any radial depen- 
dence of the field strength, so for the CCR model, the parameter 
space is defined by [D,., w, Bq,xq\. 

The best-fit values obtained by iRand & Kulkariiil ([T989h for 
the disk component were: Dr - (0.6 + 0.08) kpc , w = (3.1 ± 
0.08) kpc and Bq = (1 .6 + 0.2) pG. 

2.4. Logarithmic Spiral Arms (LSA) model 

This model was used bv lPage et alJ (l2007h to describe the distri- 
bution of the large-scale pattern of the polarization angle in the 
22 GHz WMAP 3-year data. The equation which describes the 
field is: 



where 



and 



Br - Bo smi//(r)cosx(z) 
B^ = Bo cos if/{r)cosxiz) 
Bj = Bo sin^(z) 



(8a) 
(8b) 
(8c) 



8 kpc 



Xiz) = Xotanh 



1 kpc 



Note that the LSA model is essentially an axisymmetric model 
in which the pitch angle is not constant, being magnetic field 
lines logarithmic spirals. According to our definition of pitch an- 
gle for the ASS model, iff{r) would play the role of pitch angle, 
in which we would have a constant part, i/'o, and a characteris- 
tic amplit ude for the logarithmic dependence of the arms, i^i. 
Following iPage et"ai] (|2007), Bo(r) is assumed constant, with a 
value of 3 pG. 

Thus, the parameter space is defined by [t//o,if/i,X()]. The 
proposed values fo r the different par ameters are [i/'o, (/'i,Afo] = 
[27°, 0.9°, 25°] (see'Page et al.', '20071, and also the erratum avail- 
able at the LAMBDA web sitd3). 



2.5. Bi-Toroidal (BT) model (halo model) 

Some authors have found hints of a halo or thick disc field com- 
ponent with opposite directions in both hemispheres. For exam- 
ple, iHan & Wielebinskil (l2002h and iProuza & Snifdal (l2003h de- 
tected it with maximum strengths at a large height of 1.5 kpc 
at both sides above and below t he plane, b eing the maximum 
strengths of 1 pG approximately. jSunetaTj (12008) give a more 
complete description of this double torus, being the height of 
maximum strength again at 1.5 kpc, but its maximum is much 
stronger (around 10 pG). 

Following this scenario, we propose a possible configuration 
for the components of the magnetic field which do contain a dif- 
ferent sign in both hemispheres, and it is given by: 



B, = 



B~ - constant 











|exp 


[20-2 j 



(9a) 
(9b) 
(9c) 



where cri and <T2 are two constants (measured in kpc) which 
encode the characteristic scales of variation of the field with the 



" http://lambda.gsfc.nasa.gov/product/map/curTent/map_bibIiography.cfm 
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vertical distance, and do take into account in a simplified way 
the change in the sign. For our computations, we fix the B, value 
to be 0.2 yL(G (Han & Oiao, 1994), and we only consider the case 
of a radial variation of Bo{r) given by equation (|5]l. Thus, the 
parameter space for this model is specified by [r[, 0-1,0-2]. 

3. Methodology 

In order to obtain constraints on the different parameters for each 
one of the models described above, we have performed a system- 
atic comparison of the predicted polarized intensity due to syn- 
chrotron emission with the observed map at 22 GHz by WMAP 
satellite. Here, we describe the relevant details of the dataset, the 
numerical procedure to compute the synchrotron maps for a cer- 
tain GMF model, and the model selection criterion that we have 
adopted for our analyses. 

3.1. Description of fhe K-band WMAP5 data 

The analysis of this paper is based on a comparison with the 
K-band (equivalent to 22 GHz) polarization m ap obtained by 
the W MAP satellite after five years of operation (iHinshaw et aL . 
12009 ). This map is pubUcly available in the LAMBDA website^ , 
and it is given in HEALPdO format (Gorski et al., 2(^1). 

Figure[T]shows the all-sky Stokes Q and U maps at 22 GHz, 
degrading the resolution to enhance the large-scale pattern, using 
a nside =16 pixelization (which corresponds to a pixel size of 
3.66°). For all the computations in this paper, we will use these 
degraded maps as the input data. As shown by the WMAP team 
dPage et al.L l2007h . at this frequency (22 GHz) the large-scale 
pattern observed in polarization is completely dominated by the 
galactic contribution, and the CMB component is sub-dominant. 
Thus, in our analyses we can safely neglect the contribution of 
the CMB to the polarization map. 

From these two observables (Stokes's Q and U parameters), 
one can obtain the map of the direction of the polarization angle 
(hereinafter, PA) as; 

r<„.(«) = ^arctan(^j + ^, (10) 

where n is the direction of the line of sight. Note that we in- 
clude in our definition the n/2 factor, so equation ( fTOl l represents 
the angle of the magnetic field, and takes values in the [0, n] re- 
gion. Note that the polarization convention adopted here is the 
one described in HEALPix, labelled as COSMO, which differs 
from the usual lAU convention in a minus sign for Stokes U pa- 
rameter. In addition, the WMAP definition of Stokes parameters 
includes an additional 1 /2 fac tor with respect to the definition 
used by IChandrasekhaH (1 19601) . which will be the one adopted 
here. All these quantities (U, Q and PA) are defined in a galac- 
tic (heliocentric) coordinate system. This has to be taken into 
account when comparing the observed maps with the models. 
Figure |2] shows the observed direction of the PA at the same 
3.66° (nside - 16) resolution. 

3.1.1. Noise maps 

In order to perform the model selection, we also need to estimate 
the noise maps associated to the data. It is important to note that 
those noise maps should account for the covariances introduced 

^ http://lambda.gsfc.nasa.gov/product/map/current 
http://healpix.jpl.nasa.gov/ 



Q (22GHz) 




-0.10 ^^^^H O.li) mK 



Fig. 1. Observed Stokes Q (top) and U (bottom) maps at 22 GHz 
from WMAP5 data. The maps are shown in a Mollweide projec- 
tion of HEALPix in Galactic coordinates. The center of the map 
corresponds to (/, b) = (0°, 0°), and the graticule increases with 
Al - Ab = 45°. The maps are degraded to nside = 16 to enhance 
the large scale pattern. Units are mK (thermodynamic). 



PA (22GHz) 




^^^^^^^Hb^ .^^^^H 179 deg 



Fig. 2. Observed direction of the polarization angle (PA) at 
22 GHz, obtained from the two maps shown in Fig. [1] PA is de- 
fined here as the local direction of the magnetic field (see equa- 
tionfTOli. Units are degrees. 



by all components which are present in the observed data but are 
not included in the theoretical model. In particular, it should ac- 
count for the instrumental noise component and well as for the 
random component of the GMF which is not included in our the- 
oretical model. The first one could be easily estimated based on 
the information provided by the WMAP team about instrument 
sensitivity and the overall integration time that the satellite has 
spent on each particular pixel. However, this would not account 
for the second part of the covariance. Thus, in order to model 
all the different contributions to the covariance, we follow a dif- 
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ferent procedure which makes use of the fact that the original 
WMAP-K band maps have a much better angular resolution than 
the pixel size which is finally adopted in our analyses. 



Noise maps for Q and U. In this paper, we have probed three 
different methods to characterize the noise maps associated to Q 
and U. In addition to the pure instrumental noise, all these three 
methods attempt to estimate the contribution to the total covari- 
ance matrix of the residual astrophysical components which are 
not included in our modelling (for example, the random compo- 
nent of the magnetic field, or the variance introduced by point-to- 
point variations of the mean level of the galactic emission). All 
methods produce a noise map at nside = 16 resolution, which 
provides a pixel size of ~ 3.66°. The three methods to obtain the 
ctq map (or equivalently, the ctu), are: 

- M ethod 1 . This corresp onds to the same procedure described 
in lJansson et al.l (l2008l) . We start from the observed WMAP 
Stokes-Q map at full resolution (nside = 512), degrading it 
to ~ 0.92° pixel resolution (nside = 64). For a given pixel 
i within our nside =16 pixelization scheme, we obtain the 
associated noise a-Q{i) by computing the square root of the 
variance of the ~ 0.92° pixels inside a radius of 2° from the 
center of our pixel /. 

- Method 2. We start from the observed Q-map at nside = 512, 
and we convolve it with a Gaussian of FWHM= 1°. For a 
given pixel / within our nside =16 pixelization scheme, the 
associated noise crQ{i) is computed as the square root of the 
variance in each pixel of the smoothed map within a radius 
of 2° from the center of our pixel /. 

- Method 3. We start from the observed WMAP Stokes-g map 
degraded at nside = 16. For each pixel the noise crg(i) is 
computed by obtaining the square root of the variance in- 
side a circle of radius r ~ 7.4° (i.e. twice the pixel radius). 
This method produces an unbiased estimate of the noise map 
in the case of uncorrected noise, provided that the scale of 
variation of the noise map is larger than 7.4°. 

Each one of these three methods produce different results 
for the noise map in certain areas, specially those dominated 
by galactic features. In turn, these differences imply significant 
changes in the goodness-of-fit statistic (up to a factor of 2 - 3 in 
some cases). In general, we find that the larger noise amplitudes 
are obtained with method 1, while the smaller amplitudes are 
obtained with method 2. Thus, and as a conservative approach, 
we decided to present the computations of this paper using the 
method 3. The implications of this uncertainty in the determi- 
nation of the noise maps are discussed in section |5] Noise maps 
for (Q, U) parameters, obtained with method 3, are shown in 
Figure [3] These maps have been used to obtain the final results 
presented in Tables[3]and|4] 




1.9 ^^^^^B- .^^^B SB. 4 deg 



Fig. 3. Noise maps for the WMAP5 datasets presented in Fig.[T] 
and |2] It is shown the Stokes Q (top). Stokes U (middle) and 
position angle (PA, bottom) maps. Units for Stokes's Q and U 
maps are mK. Units for the PA map are degrees. The noise map 
for PA has been computed as fluctuations around the observed 
WMAP5 K-band map. 



Noise maps for PA. Once these two noise maps {ctq and cry) 
have been obtained, one can in principle obtain the noise map 
associated to the PA (cTy) from them. However, there is an im- 
portant point to stress here. Equation (fTOl l. which defines the PA, 
is not linear in Q and U. This fact implies that the variance map 
for the PA will depend on the particular model which is used 
to compute the average value; or in other words, the noise map 
for PA will depend not only on ctq and cru, but also on Q and 
U themselves. Therefore, when doing the model selection, the 
noise map for PA will be also a function of the model. 

For illustration, we compute here the noise map for PA, us- 
ing as reference model the observed WMAP-K map. We use 



here a Monte Carlo method, drawing A^sim realizations of pairs 
of {Q, U) maps, with mean equal to the observed (Q, U) maps, 
and variance given by ctq and cru, as computed in previous step. 
We have checked that using A^sim = 5,000 noise realizations is 
enough to get convergence on the a-y map at the per cent level. 
Figure [3] also shows the noise map for the PA, obtained for this 
particular case. 
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3.2. Producing the synchrotron polarized emission for a 
given magnetic field model 

To obtain the predicted polarized synchrotron emission for any 
of the Galactic magnetic field models described above, we have 
developed a code which only includes the relevant physics, and 
it is optimised in terms of computational time with the aim of ef- 
ficiently perform all the computations. The code works directly 
within the HEALPix pixelization scheme, and obtains the pre- 
dictions for the synchrotron emission directly at the resolution 
level that we have chosen (i.e. nside - 16). As discussed above, 
this resolution is enough as long as we are interested in the large- 
scale pattern of the Galactic emission. 

In general, assuming that the cosmic ray spectrum is a power 
law distribution of spectral index p, we can predict the Stokes's 
parameters that characterize the polarization of the synchrotron 
emission at a certain frequency by computing the emissivity (en- 
ergy per unit time per unit volume per frequency per solid an- 
gle) in the two orthogonal directions, parallel and perpendicular 
to the projection of th e magnetic field on the plan e of the sky. 
Following the same as iRvbicki & Lightmanl(ll986h . we have: 

/ N N (A-nmcV^ \_p . £+1 

e±(v) = N{r,z) ^ , I ^ ' (Bsma) 2 



U 12 



Sttotc^ \ 3e 

.I'+i 



p + l 



4 12/ U 12 



(11) 



where we explicitly introduce a minus sign in the equation 
for the U component, in order to convert from the lAU con- 
vention for the polarization to the HEALPix convention which 
is used in the WMAP maps. The Ku(v) and Kq{v) constants 
also include the conversion factors between brightness and tem- 
perature. At 22 GHz, we can safely use the Rayleigh-Jeans 
approximation. Substituting the numerical values, we obtain 
Kq(v) = 1.41 X 10" mK cm^ (jiGy'- kpc K and Ku(v) = 
1.25 X 10" mK cm^ (jjGT^ kpc The simulated PA map is 
derived from these two equations ( fTSl ) as: 



y(h) = 0.5 arctan 



( -Ku(v) NeimB.Bydz 



+ 5 (16) 



where Z?v and By represent, in our coordinate system, the two 
components of the magnetic field which are perpendicular to the 
line of sight. 

Finally, when predicting the expected synchrotron emission 
for a particular model, we have included an additional restriction 
on the line-of-sight integration, by excluding those points whose 
galactocentric radial coordinate is smaller than 3 kpc or larger 
than 20 kpc. The first restriction excludes the inner region of 
the Galax y, where large devi ations from the regular pattern are 
expected (iLa Rosa et al.Ll2006,) . while the second one introduces 
a radial cut-ofi^. In any case, we have checked that the results are 
robust against changes in these numbers. 



^ VSe-' l4-nmc\ - , . ^p+i 



1 

12 



2-^P^19 
4 12 



p + l 



(12) 



where B is the magnetic field, v is the frequency, and e and m are 
the electron charge and mass, respectively. The function N(r, z) 
represents the electron number density at the corresponding po- 
sition (r,z) in the Galaxy, and it is obtained from equation ([T]i. 
From these two equations, the polarized intensity at a given fre- 
quency is obtained by integrating the emissivity along the line of 
sight: 



Iy(Z,h)^ j 



[eAv, z, h) - e||(v, z, h)] exp^'^^^^-"' dz (13) 



where we have set the coordinate system in such a way that z- 
axis represents the line-of-sight direction, and the other two di- 
rections are contained in the plane on sky, with y-axis point- 
ing east and x-axis pointing south (i.e. HEALPix coordinate 
convention, as explained in Sect. 13. 11 1. With t hese definitions, 
the S tokes's Q and U parameters are given by dChandrasekhaii 
fT%(D : 



Gv = /.v - ly 

Uv -2 -\/7^ -^cos 6 



(14a) 
(14b) 



where in our case, 6 = 0° (i.e. no desphase is assumed). 

Inserting equation (fT3] ) into equations ( fT4b . we obtain the 
simulated Q and U components along the line of sight (z-axis) 
by numerical integration: 

Qy(h) = KQ(v) f NAn)[Bl-Bl]dz (15 
Jlos 



J LOS 

Uy{h) = -Ku{v) \ NAn)2B,Bydz 
Jlos 



a) 
(15b) 



3.3. Exploration of the parameter space 

As described in Section |2] for each one of the families of GMF 
models, we have a set of parameters which define each particu- 
lar model. Given that in all cases, the dimension of the parameter 
space is small (there are, at the most, four parameters describing 
a particular model), we decided to carry out the exploration of 
the parameter space using a grid-based approach. For higher di- 
mensions, a Monte Carlo method would be more appropriate. 

For each one of the different GMF models, we have con- 
structed three different grids of models, which we label as "liter- 
ature", "blind" and "non-blind". The first one is centered around 
the average values which are found in the literature for each one 
of the different parameters. The second one spans the maximum 
range which is reasonably expected for each particular param- 
eter. Finally, the third one is built a-posteriori, once the model 
selection has been performed on the previous grid, by centering 
the new grid around the best-fit parameters for each case. Table[T] 
summarizes all the relevant parameters for each one of these 
three grids. In total, we have computed more than one million 
models (290,000 for the blind grid, 970,000 for the non-bhnd, 
and 51,000 for the literature) for all the different GMF models 
described in Section|2] Each one of those models corresponds to 
a set of three maps (Q, U and PA) of the expected synchrotron 
polarized emission of the sky at 22 GHz. 

As indication, the average execution time in a standard desk- 
top computer for the computation of a particular model requires 
^ 4 seconds of CPU time. Thus, the total CPU time for the con- 
struction of all the grids is around 1,500 CPU hours. 



3.4. Model selection and parameter estimation for each GMF 
model 

Once we have explored the parameter space with these three 
grids, we have derived the best-fit parameters for each one of 
the GMF models using a bayesian approach. To this end, we 
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Table 1. Exploration of the parameter space. For each GMF model, we show the range of values which has been used to produce 
three grids (see text for details). For each parameter, the three values indicate the minimum, the maximum, and the step size 
(uniform) which was used to build the grid. 



Model 


Parameter 


Blind exploration 


Non-blind exploration 


Literature 


LSA 


ill (°\ 
XoC) 


10, 80, 2 
10 10 OS 

iU, iU, KJ.J 

0, 40, 1 


50, 75, 0.5 
10 ions 

iU, iU, U.J 

15, 50, 0.5 


30, 40, 0.2 
IS 1 S 9 

i . J, i . J, U.z, 

15, 30, 0.5 


ASS(const) 


BoijiG) 

P{ ) 

Xo(°) 


0.5, 8, 0.5 
0, 40, 2 


0.2, 9.8, 0.2 
0, 40, 0.5 


1,3,0.5 
1 1 n ^ 

0, 20, 1 


ASS(r) 


riikpc) 
pC) 


0, 60, 2 
-30, 30, 1 
n 40 9 


0.5,60.5,1 
0, 30, 0.5 
so s 

U, JU, U.J 


0, 20, 1 
-15, 15, 1 
90 1 

U, Z.U, 1 


BSS±(const) 


pn 

/10(, ) 


0.5, 8, 0.5 
-30, 30, 1 
40 2 


0.5, 8.5, 0.5 
0, 40, 0.5 
^S s 

U, JJ, U.J 


1,3,0.5 
-15, 15, 0.5 
40 2 


BSS±(r) 


riikpc) 
PC) 
XoC) 


0, 60, 2 
-30, 30, 1 
0, 40, 2 


0.0, 60.0, 1 
0, 30, 0.5 
0, 35, 0.5 


0, 20, 1 
-15, 15, 1 
0, 20, 1 


CCR 


Drikpc) 
w(kpc) 
BoQjG) 
XoC) 


0.1, 10.1, 1 
0.1, 19.1, 1 
0,9,1 
0, 40, 2 


1, 11,0.5 
3.1, 20.1, 1 

2, 10, 0.5 
4, 50, 1 


0,1.5,0.1 
2,4,0.1 
1,3,0.2 

10, 40, 1 


BT 


riikpc) 
CTiikpc) 
o-iikpc) 


0.5, 60.5, 2 
0.01, 10.01,0.5 
0.01, 10.01,0.5 


0, 50, 1 
0.01,5.01,0.02 
0.01,20.01,0.5 





tions are obtained, and then marginalisecQ over all the unwanted 
parameters. At the end, we end up with marginal probability dis- 
tribution functions for each one of the parameters. From these, 
confidence intervals are derived as the 0.16, 0.5 and 0.84 points 
of the cumulative probability distribution function. Thus, our 
parameter estimate is the median of the marginalised posterior 
probability distribution function, and the confidence interval en- 
compasses 68 per cent of the probability. 

3.4.2. Priors 

For the analyses in this paper, we have not introduced prior in- 
formation on the parameter values describing any of the mod- 
els. This is equivalent to say that we have implicitly adopted a 
top-hat prior in all the parameters, where the top-hat function 
is defined by the parameter ranges presented in Table [1] Thus, 
in all cases, the evaluation of the posterior would reduce to the 
computation of the likelihood function (X.)- 

However, and for the case of (Q, U) analysis, we have 
slightly modified the standard analysis in the following way. As 
discussed in Section |2] the amplitude of the CR electron spec- 
trum in the solar neighbourhood is highly uncertain. This will in 
turn imply a large uncertainty in the recovered strength for the 
magnetic field, and moreover, it might produce a bias on the re- 
covered GMF amplitude. To account for this additional degree 
of uncertainty (at least at first order), we have introduced an ad- 
ditional parameter e, which multiplies to both the predicted Q 
and U maps for a given GMF model. Note that such parameter 
have no impact on PA. If the CR electron density were entirely 
correct, then we would have € - 1 . If there is an uncertainty in 



The marginal distribution functions are obtained by integrating the 
joint distribution function over the variables being discarded. For exam- 
ple, in the case of the LSA model, the marginal distribution function for 
-£(iAo) is derived by integrating the joint distribution L(4/o, i^i ,X(i) over 
ipi andxo- 



have to both compute the likelihood function {£.), and to pro- 
vide an expression for the priors. Once we have obtained the 
best-fit parameters for each GMF model, different models with 
be compared in terms of the reduced -statistic. 

3.4.1. Likelihood function 

In general, we may assume that the likelihood function is de- 
fined as a multivariate Gaussian when it is written in terms of 
the observables, i.e.: 

lnX = -i;r'- (17) 

If we assume that the coiTelations between the different pixels 
are negligible, then we have: 

;^' = Z^ (>") 



where x,- represents the observational data, k, the simulated data 
and cr^, the associated noise covariance. In our case, we have 
performed two different evaluations of the likelihood. 

- The first one coiTesponds to a direct comparison of Stokes's 
(C U) parameters. In this case, we have / = 1, 2Npix, and 
Xi = Qi for / = l,Npi/, and Xi = t/;-Wp„ for i = Npi^ + 
1, ...,2Npix- This case will be noted &^Xq u- 

- The second case corresponds to the comparison of PA, so 
now we have / - \,...,Npix, and jc, = y,. Note that in this 
case, cTj will depend on hi, and thus the PA noise map needs 
to be computed for each particular model. This case will be 
noted as Xp^i and will be used for comparison with the re- 
sults of the previous case. 

Once these two functions (y^ and x]>j<) are evaluated in 
all the data-points of the different grids, the posterior distribu- 
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this parameter due to the modelling of the CR distribution, we 
account for it by introducing a Gaussian prior for this additional 
parameter; 

-2ln£Qu=^-^^+xiu (19) 

and we marginalise over it. Given the existing uncertainty on 
the CR electron density, we have chosen a conservative value of 
CTf - 0.8. The marginalisation over this additional parameter, e, 
can be done analytically, yielding: 



In^Qf/ = 



+ 2C InB 
4-B ~ 



where: 



1 

2^ 



C = 



1 

ZXiki 1 
. ~^2^ 



2 ^ 2cri 



(20) 



(21) 



(22) 



(23) 



Note that this scheme is completely equivalent to the marginal- 
isation over calibration un certainties which i s also adopted in 
CMB experiments (see e.g. lBridle et al.Ll2002l) . 



Table 2. Galactic masks used in the analyses. Columns 1, 2 and 
3 provide the mask identification number and the mask descrip- 
tion. Column 4 shows the total number of available pixels for 
the Q-U analysis. Note that for the case of a PA analysis, we 
would have 1 /2 of this value. Last column shows the sky fraction 
available after applying each mask. 



ID 


Mask definition 


Region probed 


N,„, (QU) 






(regions excluded) 






(%) 


1 


GC 


global 


5432 


88.4 


2 


GC + Loops + mask WMAP 


global 


3452 


56.2 


3 


GC + Disk 


halo 


4416 


71.9 


4 


GC + Disk + Loops + mask WMAP 


halo 


2300 


37.4 


5 


GC + Halo 


disk 


1016 


16.5 


6 


GC + Halo + Loops 


disk 


482 


7.8 




3.5. Masks 

Nearby structures in the Galaxy (e.g. supernova remnants) might 
distort the regular pattern of the GMF as seen in the synchrotron 
polarized emission, thus biasing the determination of some pa- 
rameters of the GMF model. To probe the robustness of the es- 
timates as a function of the region which is used for the fit, we 
have considered different masks to constrain the global, halo and 
disk components of the GMF. In all cases, we exclude the galac- 
tic center region defined as the region in which the line of sight 
crosses a region with radius 3 kpc in cylindrical coordinates, that 
it is not taken into account in our integration scheme. To con- 
strain the global component we considered two masks, listed in 
decreasing order of covered sky fraction: 

- Mask 1 . It corresponds to the Galactic center region. 

- Mask 2. It combines mask 1 with the well-known local re- 
gions with a strong polarized intensity that could be inter- 
preted as part of the polarized intensity produced by the reg- 
ular GMF . Here we have c onsidered the four "loops" de- 
scribed in lBerkhuiisen et aP (T971). and we have added the 
polarization mask from WMAP team, which includes a bet- 
ter masking of the North Galactic spur region and several 
other small objects (e.g. LMC). 

To constrain the halo component we exclude the following 
regions: 

- Mask 3. It excludes the disk defined as the emission contains 
in \b\ < 10°. 

- Mask 4. It excludes the emission of the disk (mask 3) and the 
"loops" defined in mask 2. 

For completeness, we have also considered two masks to con- 
strain the field pattern in the disk. There are: 

- Mask 5. It excludes the halo region defined as the emission 
obtained for 1^71 > 10°. 



Fig. 4. Regions used for the definition of the six masks adopted 
for the analyses (see text for details). 

- Mask 6. It excludes the halo and the "loops" described in 
iBerkhuiisen et al. l (fT97l . 

In Figure |4] we show all regions which have been described 
in this subsection. Figure |5] shows the polarization mask used 
by the WMAP teanfl Table |2] presents the detailed information 
about the sky coverage of each one of these masks. It is also in- 
dicated the number of available pixels for the analysis (i.e. num- 
ber of terms in the summation in equation ( fTsT l. This quantity is 
relevant in order to compute the reduced ;i'2 for the best-fit mod- 
els. As a reference, note that in this pixelization (nside - 16), a 
whole-sky map contains 3072 pixels. Therefore, we would have 
in this case NpiAQU) - 6144. 

We would like to mention that, even if we use these masks 
to eliminate the effects of the random local spurs, other unde- 
tected random local features could be a source of errors in our 
large-scale models and these errors are difficult to quantify. But 
probably our spur masks eliminate the major contribution of lo- 
cal features. 



4. Results and discussion 

For each one of the masks described in Table|2] and for each one 
of the GMF models presented in Section|2] we evaluate the pos- 
terior distribution in each one of the three grids, and marginal- 
ising over the relevant parameters, we obtain the corresponding 
confidence regions. Our default analysis uses the Q-U maps, and 
the noise maps obtained with method 3 in section [J. 1.1 1 The re- 
sults are summarized in Tables [3] and|4] 

In order to evaluate which model better reproduces the data, 
we have used as a goodness-of-fit the reduced ;i'2 statistic, which 



http://lambda.gsfc.nasa.gov/product/map/dr3/masks_get.cfm 



10 



Ruiz-Granados et al.: GMF and WMAP5 data 




Fig. 5. Mask for the galactic polarized emission used by the 
WMAP team. 

is obtained as the ratio of the the minimum value for x^(= 
-2 In -Cqu) to the number of degrees of freedom (hereafter dof). 
The number of dof is obtained as A^pu - M, being Npix the num- 
ber of terms in Equation ( fTSl ). and M the number of parameters 
for the considered GMF model. Last two columns in each table 
present the values for the minimum and the reduced 

4.1. The magnetic field in tiie Galactic halo 

The results for the halo field are summarized in Table [3] 

Our reference computation corresponds to the case labelled 
as mask 4, in which we exclude the emission of the disk, 
the "loops" (including the polarization mask provided by the 
WMAP team), and the galactic center region (which is not ac- 
counted for in our analysis). For this reference case, the model 
having the minimum reduced is the ASS(r), although the 
other two axisymmetric models (LSA and ASS) provide approx- 
imately the same goodness-of-fit. Thus, the large-scale 22 GHz 
polarized synchrotron emission seems to be more compatible 
wi th some type of axisymmetry, a conclusion also reached 
bv lPage et al.l(|2007[) . and also compatible with the results shown 
by ISun et al.l (12008*). For illustration. Figure |6] shows the field 
pattern of the best-fit ASS(r) model at z = 4 kpc; figureQshows 
the marginalised one-dimensional posteriors distributions for the 
parameters of this model; and figure[8]shows the predicted Q, U 
and PA maps for the same best-fit model. We now discuss sepa- 
rately each one of the relevant parameters. 

Radial scale. For this ASS(r) model, the derived constraint on 
ri is < 2.5 kpc (95% confidence level). This parameter essen- 
tially controls the distance at which the magnetic field is no 
longer constant and begins to decrease proportional to r ' (see 
Eq. |5]l. The obtained value is indeed very small compared with 
the radial scale of the electron density or any other scale distance 
in our galaxy, suggesting that the data require an important vari- 
ation of the field in the inner part of the halo, probably due to the 
pres ence of stronger magnetic fields at the galactic center (see 
e.g. IRov et a\!. '2008*). Indeed, in the literature, the halo model 
proposed in Prouza & Smida (20Q2) requires also a small radial 
scale of ri =4 kpc, being the radial dependence oc j- exp (--^Tp) 
in that case. 

Field strength. For the ASS(r), we adopted a fixed values for 
the magnetic field strength of Bq = 3 fiG at th e solar neighbour- 
hood. The original LSA model proposed by iPage et al.l (l2007h 



also assumed this fixed value. However, it considered as a free 
parameter for the ASS model, and was found to be in agreement 
with that value (Bq = 2.8+2-^ /iG). 

Pitch angle. The pitch angle was considered in the ASS(r) and 
ASS models as a constant free parameter In both cases, we con- 
sistently obtain a value between 24° and 26° respectively. For the 
LSA model, it is considered as a radial function with a logarith- 
mic dependence, and in this case the derived pitch angle at the 
solar neighbourho od is again 26°, a value which is consistent 
with that given bv iPage et al.l (l2007h of p — 27°. As discussed 
below, these values are larger than the typical pitch angles ob- 
tained for the disk field. 

Tilt angle. This parameter controls the vertical structure of the 
field. The derived tilt angle in all the axysimmetric models is 
of the order of 30°, which implies a vertical component close 
to 1 juG at z = 1 kpc. This vertical field component could be 
identified as the poloidal component corresponding to the dipole 
fiel d responsible of the mG vertical component at the very cen- 
ter jHan & Oiaoll 19941) . It could correspond as well to the ver- 
tical component of the cluster field diffused into the disk by tur- 
bulent magnetic diffusion (e.g. Battaner & Florido.,2000,) . in this 
case to the Local Group field. 

Other models. Table [3] show a good consistency between the 
three axisymmetric models. The rest of the families considered 
here provide slightly worse reduced figures, although any of 
them can not be clearly rejected. In general, the two families 
of bisymmetric models provide a slightly poorer goodness-of- 
fit. The derived radial scales, field strengths and pitch angles are 
in general similar to those found for the axisymmetric models. 
However, the tilt values (xo) are considerably lower, and in some 
cases negligible. This could be due to a compensation produced 
by the inherent reversals of the fields for a given direction in a 
bisymmetric configuration. 

The two remaining models (CCR and bi-toroidal) produce 
the poorer fits (x^ - 2.1), but again they can not be rejected. 
Attending to this values, the CCR points to the existence of a 
reversal at r ~ 3 kpc from the galactic center. The bi-toroidal 
model has been suggested as a possible explanation for the halo 
dou ble torus th at it is interpreted by iHan et alj d 19971 1 19991) 
and IHmI (120091) as a consequence of an or - Q effect. Another 
possibility is that the aforementioned vertical field diffused from 
the galaxy cluster could be twisted by differential rotation in 
the vertical direction, producing toroidal fields above and be- 
low the plane with opposi te directions in both hemispheres (e.g. 
iBattaner & Floridoll2000l) . 

Effect of the loops and the disk emission on the determination 
of the halo field As described in section [331 nearby structures 
in the Galaxy might introduce biases on the recovered parame- 
ters of a given GMF model. Moreover, the emission of the disk 
could contaminate the halo field as discussed in section [1] If the 
disk emission is not excluded by the mask, we could constrain 
the global component of the GMF and quantify the impact of 
the disk emission on the parameters describing halo field. The 
corresponding results are also shown in Table [3] 

The influence of disk emission can be seen by comparing 
the results from masks 2 and 4. In both cases, the best-fits are 
obtained for axisymmetric models but when disk emission is not 
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Table 3. Best-fit parameters for the global and halo component of the GMF model (masks 1, 2, 3 and 4), based on the analysis of the 
Stokes's Q and JJ parameters. For each mask (see TablelJ]) and each GMF model (see section|2]i, it is shown: (a) confidence intervals 
for the parameters listed below derived from the cumulative probability distribution function, obtained after marginalising over the 
rest of parameters; (b) minimum;t'^ = -2 In Hqjj for the "best-fit model", defined as the one having the smaller and (c) minimum 
divided by the number of degrees of freedom. Confidence intervals encompass the 68% of the probability, except in those cases 
in which only upper (lower) limits can be given, in which we provide the 95% limit. For a quick reference, the parameters which 
define each GMF model are: LSA (i/'o,(Ai,;ro); CCR model (D„ w, Bo,;t-o); ASS (Bo,p,A'o); BSS± (Bo,p,A-o); ASS(r) (ri,/?,;to); 
BSS±(r) (ri , p,x<i)\ BT {r\ ,(T\,(T'i). Our reference computation corresponds to the one labelled as mask 4. 



Mask Model Confidence regions MinO^'') for best fit Min(:^^)/d.o.f. 

i LSA 23.9+;;-5, > 5.8, 32.4 ± 1.1 10668.6 L9 

1 CCR 8.5±0.2r<5.1,< 8.6, 20.6!JJ-^ 13828.8 2.5 

1 ASS 2.8+J' 23.9!5!-^, 20.5+U-^ I069I.I 1.9 

1 BSS+ 2.8!^;^, 24.5+f ^, < 2.0 12056.9 2.2 

1 BSS_ 2.9!jf, 25.3 ± 0.4, < 2.0 11964.6 2.2 

1 ASS(r) > 38.0 , 23.8+]; 5, 21.7 ±0.7 10703.0 1.9 

1 BSS+(r) > 50.8, 24.4 ±0.6, < 2.0 12103.4 2.2 

1 BSS_(r) >42.4, 25.1 ±0.4, < 2.0 11984.8 2.2 
\ BT > 52.5, < 1.0, > 1.1 15048.3 2.8 

2 ESA 23.7 ±0.5, > 8.1, 22.0+|;-^ 54563 LS 
2 CCR 2.5±0.2, <5.1,<8.6, 13.r±0.9 7053.3 2.0 
2 ASS 2.8!^-^, 23.3!jJ*, 22.4 ± 0.9 5504.5 1.6 
2 BSS+ 2.7!;^, 24.6+°' < 2.0 6717.1 1.9 
2 BSS- 2.7+J-^,24.5+i|-^,< 2.0 6793.8 1.9 
2 ASS(r) > 42.5,23.1 ±0.5, 23. 8+° 5 5523.6 1.6 
2 BSS+(r) > 47.2, 24.8+];-^, < 2.0 6751.0 2.2 

2 BSS_(r) > 43.2, 24.6+J5, < 2.0 6820.0 1.9 
2 BT >42.0, < 1.0. 0.48 ±0.2 7641.4 2.2 

3 LSA 25.5 ± 0.4, < -4.9, 25.9 ±0.8 87T48 O 
3 CCR 8.5±0.2, <5.0, < 8.6, 22.2;| 'J 11735.6 2.6 
3 ASS 2.8+J", 25.4.0 ± 0.4, 25.4!jj| 8729.8 1.9 
3 BSS+ 2.7+i^, 24.6 ± 0.4, < 2.0 9741.2 2.2 
3 BSS_ 2.8+1-^, 25.4+"-^, < 2.0 9623.0 2.1 
3 ASS(r) 2.3 ± 0.6, 23.7+[j|, 29.4+j;* 8533.2 1.9 
3 BSS+(r) 20.5+^'„'', 24.2*5, < 2.0 mS.\ 2.1 

3 BSS_(r) < 2.2, 24.9+°*^^ < 2.0 9478.4 2.1 
3 BT > 33.8. 2.9+j;-°, > 4.7 12099.3 2.7 

4 LSA 26.0+^-^, > -6.7, i2A*\\ 3mO L4 
4 CCR 3.0 ± 0.2, 4.0 ± 0.4, < 8.7. 20'.6+] f 4782.9 2.1 
4 ASS 2.8+JJ. 26.0+°-^, 32.5+|-? 3338.6 1.4 
4 BSS+ 2.6;i*, 23.4 ±0.7, < 2.0 4164.6 1.8 
4 BSS_ 2.5+1*, 23.4+];-5, < 2.0 4153 \ % 
4 ASS(r) < 2.5, 24.3 ± 0.6, 30.3+J i 3195.7 1.3 
4 BSS+(r) < 2.0, 24.0 ± 0.6. < 2.0 4037.1 1.7 
4 BSS_(r) < 2.0 , 23.8 ± 0.6, < 0.2 4017.1 1.7 
4 BT > 33.8, 2.5 ±0.2. > 3.6 4824.9 2.1 



masked out, the reduced-;i^^ becomes slightly poorer In general, 
the pitch and tilt angles and the field strength remain unchanged, 
with the important exception of the radial scale factor for ASS(r) 
and BSS±(r) models. The inclusion of the disk emission in the 
analysis increases drastically this radial scale, probably because 
as shown below, the disk does not require strong radial varia- 
tions. 

Finally, we can also evaluate the impact of the loops on the 
fit by comparing the results of masks 3 and 4, or 1 and 2. The 
basic conclusion in this case is that including the loops regions 
in the analysis do not bias significantly the results (even for the 
radial scale parameter), but the quality of the fits get worse in all 
cases. 

4.2. The magnetic field in the disk 

For completeness, in this work we have also used two masks 
(5 and 6) to study the magnetic field in the Galactic Plane, by 
masking the halo emission. Table |4] summarizes the constraints 
on the different parameters for this case. The reference mask 
now is number 6, which also excludes the contribution of loops. 
However, the available sky area for the fit in this case is very 
small (7.8% in total, which corresponds to approximately a 40% 



of the total area of the disk). Because of this limited area, the 
conclusions on the magnetic field parameters might be uncer- 
tain. Nevertheless, we consider that still it is important to com- 
pare these results with the numbers obtained with other methods. 

Focusing on the mask 6 alone, the three axisymmetric mod- 
els (LSA, ASS and ASS(r)) provide practically identical value of 
the goodness-of-fit, which is slightly better than the other cases. 
For these three models, the derived pitch angle values are lower 
than in the halo case. In the solar neighbourhood, the pitch angle 
of the spiral arm is ~ 18° for the stars and ~ 13 ° for all gaseous 
component^] (see IValleel 119951: IValled |2002|) . Therefore, our 
best-fit suggests that the magnetic arms fo llow the gas st r ucture 
(p ~ 14° - 15°). However, we note that iJansson et al.1 (l2009h 
fo und a valu e of p ^ 35° for the ASS+RING model proposed 
bv I Sun etal ] (I2008h . which is not compatible with the one ob- 
tained here. 

Among the two families of bisymmetric models, both of 
them provide good results, and again, the pitch angle values 
are low. Moreover, the derived constraints (between 7° - 10°) 
are fully compatible with the results obtained by other au- 
thors. In the literature, values also range from 7.2° to 11° 

' Note that these values are translated into our sign convention for 
the pitch angle. 
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Table 4. The same showed in Table|3]for disk component of the GMF model, from the analysis of the Stokes's Q and U parameters. 
Our reference computation corresponds to the one labelled as mask 6. 

Mask Model Confidence regions Mm(x~) for best fit Min(;^-)/d.o.f. 



5 ESA 16.9;; ^, -1.6+5 J, < 15.9 1729!2 TT 

5 CCR 3.0±0.8,> 14.6, < 8.8, < 6.6 1970.4 1.9 

5 ASS S.l^J-J, 17.7 ±0.8, 5.9!] i 1664.0 1.6 

5 BSS+ 4.1+1^, 12.2 ±0.4, 18.0;i-^ 1626.1 1.6 

5 BSS_ 3.9+f g, 13.6 ±0.4, < 2.0 1775.5 1.7 

5 ASS(r) >49.1,16.8±0.8,6.1± 1.2 1696.7 1.6 

5 BSS+(r) >46.0, 11.9±0.4, 19.0!|5 1^47.0 1.6 
5 BSS-(r) > 49.2, 13.0 ±0.5, < 2.3 1812.8 L8^ 

6 ESA < 19.8, > 0.0, U.y\* 6723 TJ 
6 CCR 4.5:J-^, 11.4:2|, < 8.7,21.7^^-5 757 3 j g 
6 ASS 3.0+^;°, 15.8+1 ^_ 17.3 ± 1.5 576.8 1.4 
6 BSS+ 3.9!^°, 10.8 ±0.5, 16.4 ±2.0 751.2 1.6 
6 BSS_ 3.7+2;;, 7.8!2*, 18.0+15 774.2 1.6 
6 ASS(r) > 26.6, 14.6;i 3, 18.5+j| 680.7 1.4 
6 BSS+(r) > 36.7, 10.5 ± 0.6, 17.6+1; 761.2 1.6 
6 BSS_(r) > 33.9, 7.2+°-^, 19.2+^-f 786.5 L6^ 



20 rr 




-20 ^ , , , , - 

-20 -10 10 20 

X-Galactocentric (kpc) 

Fig. 6. Large-scale pattern of the ASS(r) model at z = 4 kpc. 
This model provides the best-fit for the halo field. 



(see iHan & Oiaol (1 19941) and EiB (1200 lb for values obtained 
with Faraday rotation of pulsars; and iHeilesI (1199 6), for val- 
ues obtained with polarized starlight). For comparison purposes. 
Figure |9] shows the pattern of our best-fit BSS+ model in the 
galactic disk (z = 0), which is similar in shape to the results ob- 
tained by other authors (see e.g. figure 5 in Han & Oiao ( 1994)). 

We note that the values derived for the tilt angle in all models 
are of the order of ^ 19°, which again imply a z-dependence of 
the field strength w ithin the disk which is compatible with those 
values observed bv lHan & Oiaol d 19941) and up to 0.4//G for the 
thin disk. 

Finally, we would like to mention that we do not expect 
changes in these results if we include a more refined treatment 
of the turbulent magnetic field in the analysis, since the polar- 
ized synchrotron emission comes from the regular patter n of the 
GMF which is located in the inter-arms regions dBeckL l2007h . 
Indeed, we find that the results obtained with the other two meth- 
ods for the noise determination (see Sect. 13.1.11 ) are fully con- 
sistent with those presented here for all masks and all parame- 
ters. We only found a small dependence of the constrained value 
for xo for axisymmetric and bisymmetric models with the three 
noise maps, but which is of the order of :$ 15 per cent. 



4.3. Comparison of the result with the PA analysis 

For comparison purposes, in this paper we also have done the 
analyses using the PA information alone. As discussed else- 
where, such analysis provides a limited amount of information, 
due to the fact that the PA is not sensitive to some parameters, as 
a constant field strength. 

In general, the results using the PA provide compatible re- 
sults to those of the QU analysis, although there is a larger num- 
ber of unconstrained parameters and a poorer goodness-of-fit. As 
illustration, the ASS(r) analysis with mask 4 gives ri < 45.9 kpc, 
p ~ 20.0°, and xo ~ 21.0° Or' ~ 2.3); while for the ASS 
case it gives p ~ 20.5°, xo ~ 19.0° and no constraint on Bq 
(X^ ~ 2.3). The pitch and tilt angles are compatible with those 
derived from the QU analysis. The best fit is given by the LSA 
with lAo > 22.3°, lAi > 8.6° and xo < 18.1° (x^ ~ 2.3). 

5. Conclusions 

In this work, we have constrained the regular galactic magnetic 
field by using the polarized emission at 22 GHz. To this aim, 
we have considered nine models of the galactic magnetic field, 
each defined by three or four free parameters and for six different 
masks to interpret the polarized maps from WMAP5. The com- 
bination of models, free parameters, and masks produce a very 
large number of simulated maps to be compared with the ob- 
servational ones in Stokes's Q and U parameters, which in turn 
provide valuable constraints to determine the three dimensional 
configuration of the magnetic field of our galaxy. 

The family of GMF models that better describes the halo 
emission is the axisymmetric one, although any of the other con- 
sidered models can be rejected based on their goodness-of-fit. 
The magnetic spiral arms have a pitch angle of p x 24°, and a 
tilt angle of xo ~ 30°, implying a strong vertical field of ~ 1 
fiG at z = 1 kpc. When a radial variation is fitted, the models 
generally require a fast variation in the inner part of the galaxy 
(ri < 2.5 kpc). 

We would like to stress that an accurate determination of the 
covariance matrix which accounts for both for the noise and the 
residual astrophysical components is very complicated. In this 
work, we have explored in detail three different methods, and 
found that they may lead to differences of a factor of ~ 2 - 
3 in the goodness-of-fit values, while the values of the best-fit 
parameters for each model do not vary significantly. In practice, 
this means that rejecting a model based on the goodness-of-fit 
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Fig. 7. One dimensional (marginalised) posterior distribution 
functions for the parameters of ASS(r) halo model (top: ri, 
middle: p and bottom:;t'()) when mask 4 is considered in the QU 
analysis. 



could be inappropriate, although the relative differences of the 
statistic between models can be used for a comparison and 
for selecting preferred models. 

We have also tried to constrain the disk field parameters us- 
ing the polarized synchrotron emission, despite of the fact that 
the fitted region is very small (~ 8% of the sky). Here, the ob- 
tained results a remarkably consistent with those obtained with 
other methods. In this case, all the considered models give a very 
similar goodness-of-fit, with a (very small) preference for an ax- 
isymmetric model. The data does not require a radial dependence 
of the strength, and the pitch angle values are much smaller than 
in the halo case. Indeed, the constrained pitch angles are com- 
patible with those found by using other observational methods as 
Faraday rotation of pulsars (iHan et al.L l2006h However, we note 
that this conclusion has not been reached by recen t results where 
polar ized synchrotron emission is used (see e.g. iJansson et al.L 
120091) . The tilt angle in the disk is^o ~ 17°, which implies a ver- 
tical field structure being ~ 0. 1 iiG atz = 200 pc. This value 
is compatible with that found by iHan & Oiaol d 19941) . based on 
the rotation measure of pulsars. 




Fig. 8. Best-fit for the halo field by using Q,U parameters and 
PA excluding the disk, the galactic center and the "loops" (ASS 
model with radial dependence). 



We remark that there are still some important uncertainties in 
the modelling of the synchrotron emission. In particular, proba- 
bly we need a better knowledge of the distribution of cosmic 
rays, which is today a clear source of indeterminations in this 
type of analyses. The detailed modelling of the halo field could 
have an influence on the cosmic rays trajectories and could be 
crucial for the direct detection for the primordial magnetic fields. 

Finally, we expect that the higher sensitivity and angu- 
lar res olution in the polarized c hannels o f the PLANCK tele- 
scope dThe Planck Collaborationi 2006; Ta uber et al.Ll2010l) . flie 
low frequency channels of the QUIJOTE-CMB experiment 
(Rubiiio-Martm et al., 2010) and experiments as LOFAR and 
SKA (Beck, 2009), will provide a much larger improvement of 
our knowledge of the galactic magnetic field. 
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Fig. 9. Large-scale pattern of the best-fit BSS+ model for the disk 
emission in the galactic disk (z = 0). 
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